Depinning transition of dislocation assemblies: pileup and low-angle grain boundary 
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We investigate tlie depinning transition occurring in dislocation assemblies. In particular, we 
consider the cases of regularly spaced pileups and low angle grain boundaries interacting with a 
disordered stress landscape provided by solute atoms, or by other immobile dislocations present in 
non-active slip systems. Using linear elasticity, we compute the stress originated by small deforma- 
tions of these assemblies and the corresponding energy cost in two and three dimensions. Contrary 
to the case of isolated dislocation lines, which are usually approximated as elastic strings with an 
effective line tension, the deformations of a dislocation assembly cannot be described by local elastic 
interactions with a constant tension or stiffness. A nonlocal elastic kernel results as a consequence 
of long range interactions between dislocations. In light of this result, we revise statistical depinning 
theories and find novel results for Zener pinning in grain growth. Finally, we discuss the scaling 
properties of the dynamics of dislocation assemblies and compare theoretical results with numerical 
simulations. 
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I. INTRODUCTION 

The depinning transition of individual dislocations 
gliding on their slip plane has been widely investigated in 
the past |lcj2ip3,iA, 3\ in order to explain solid solution 
hardening |3jIE@1j that is, the increase of the yield stress 
value when solute atoms are present in a crystal. The 
presence of solute atoms changes the local properties of 
the host material, resulting in a pinning force on nearby 
dislocations d, |^ . This is not the only source of pinning, 
which can also be provided by particle inclusions or by 
immobile dislocations in other inactive slip systems [3. 
Several approximate calculations have been performed 
in the past to obtain the depinning stress from a statis- 
tical summation of individual pinning forces. In this re- 
spect, collective pinning theories have been v ery success- 
ful in the case of diffuse weak pinning forces LlhilhJ, 
whereas the theory introduced by Friedel Tsfis appro- 
priate in the case of localized strong pinning centers. 

/^From the purely theoretical point of view, disloca- 
tions provide a concrete example of a more general prob- 
lem: that of driven elastic manifolds in quenched random 
media Apart from dislocations, other examples of 

this general problem are domain walls in ferromagnets 
[15L 01 , flux hues in type II superconductors 17, 18], 
contact lines 0, , and crack fronts 0, . In re- 
cent years, a vast theoretical effort has been devoted to 



understand the depinnin g tr ansition as a non-equilibrium 
critical phenomenon [Hl^lllllllllllllllll. The mor- 
phology of a depinning manifold is generally found to be 
self-affine and can be characterized by a roughness ex- 
ponent. Other scaling exponents have been introduced 
to characterize the behavior of correlation lengths and 
times, the velocity above depinning, and the avalanching 
motion observed as the critical threshold is approached. 
Quantitative predictions of the critical exponents have 
been obtained analytic ally by the renormalization group 

^ and have been confirmed by 
numerical simulations [HlHllilllllllllllllli. In 
the course of time, a deeper level of description and un- 
derstanding of this phenomenon has been achieved, going 
far beyond a mere estimate of the depinning force, which 
has typically been the original motivation to address the 
problem. 

The analysis of the depinning transition in dislocation 
theory has often been made in the approximation of a 
dislocation line tension QjSSSIS- Hence, a disloca- 
tion is considered to be a flexible string with a constant 
line tension. This analogy is not fully accurate. In fact, 
the bending of a dislocation produces long-range stress 
and strain fields [sgIIs^ such that the self-energy of a dis- 
location line segment depends on its interaction with all 
the other segments on the line, and therefore depends of 
the particular dislocation configuration. Long range in- 
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FIG. 1: Transmission electron micrograph taken from a Cu- 
14.4at%Al single crystal deformed at room temperature; the 
image shows large regularly spaced dislocation pile-ups. Cour- 
tesy of J. Plessing and H. Neuhauser [in- 



teractions lead to a wavevector-dependent effective line 
tension. This is also well-known in the case of vortex lines 
in high temperature superconductors |38j| or, for instance, 
for a dislocation line in a vortex crystal [33 ■ In spite of 
this wavevector dependence, the main features of the de- 
pinning transition are basically unchanged , due to the 
fact that this effective line tension is only varying loga- 
rithmically with the wavevector. Numerical simulations 
indicate, however, a slight difference in the roughness ex- 
ponent in the two cases which, at present, is still not well 
understood similar surface tension approxima- 

tion is also used in the theory of Zener pinning [4l|4ll,to 
describe the interactions between grain boundaries and 
solute atoms during grain growth |42j |. 

While the behavior of an isolated dislocation pushed 
through a random distribution of obstacles is at present 
quite well understood, the results do not necessarily 
carry over to the more realistic case of collective dis- 
location motion. Dislocations interact via their long- 
range stress fields, which may induce intriguing jamming 
and avalanche-like phenomena even in the absence of 
quenched pinning centers ,4^. In most cases, one can 
not simply neglect interactions and treat dislocations as 
isolated objects. The depinning transition of interacting 
dislocation lines and/or loops of generic orientations and 
Burgers vectors in a random solute distribution seems to 
be a formidable task for theoretical treatment. In several 
instances, however, dislocations are arranged into regu- 
lar structures that are amenable to a theoretical treat- 
ment. In particular, here we analyze the depinning tran- 
sition of regularly spaced pileups and low angle grain 
boundaries (LAGB). These relatively simple structures 
are sometimes observed experimentally (see Fig. ^ and 
provide a nice illustration of the effect of dislocation in- 
teractions on the depinning transition. An early analysis 
of the depinning of a dislocation pileup was presented in 
Ref. considering explicitly the emission of disloca- 
tions from a source. 



Here we address the problem by first computing the 
stress and elastic energy associated to a small deforma- 
tion of the dislocation structure. This allows us to assess 
the validity of a local tension approximation, which turns 
out to be completely inadequate in this context. Long- 
range interactions between the dislocations involved give 
rise to much larger self-stresses than for isolated disloca- 
tions. As a result, pileups and low angle grain bound- 
aries are much stiffer than single dislocations. The elas- 
tic energy is then used in the framework of statisti- 
cal pinning theories to estimate the depinning stress. 
In particular we consider collective pinning theory and 
Friedel statistics. As an application of these results, 
we revise the theory of Zener pinning in grain growth 
which originally made use of a surface tension approxi- 
mation for the elastic self-stress. The correct treatment 
of long-range stresses and dislocation interactions inside 
the grain boundary yields a different result for the de- 
pendence of grain size on material parameters. 

At stresses close to the depinning stress, the dynamics 
of the dislocation structures we are investigating exhibits 
critical behavior which can be characterized in terms of 
scaling exponents. Using previous renormalization group 
results, we gain a complete quantitative picture of the de- 
pinning transition. In the elastic approximation, pileups 
and low angle grain boundaries are equivalent to a stan- 
dard interface depinning problem with long-range elas- 
ticity. In two dimensions, the problem can be mapped to 
a contact line or to a planar crack, which have been ex- 
tensively studied in the literature. In three dimensions, 
the self-stress is similar to the dipolar force in magnetic 
domain walls and leads to logarithmically rough deforma- 
tions. In more technical terms, d = 3 is the upper critical 
dimension for the transition, which is well described by 
mean-field exponents, up to logarithmic corrections. 

The scaling exponents associated with the depinning 
transition describe not only the morphology of the dis- 
location assembly but also its dynamics. In order to 
confirm the validity of the elastic calculations, on which 
we based the mapping with elastic manifolds, we per- 
form a series of numerical simulations for a dislocation 
pileup. We consider a two dimensional system, neglect- 
ing the deformation of single dislocations, which amounts 
to an effective one dimensional particle model. Simula- 
tions of the model display results in agreement with the 
theory and allow to illustrate some interesting dynamical 
effects. In particular, the pileup displays a zero tempera- 
ture power law creep relaxation which can be interpreted 
by scaling relations. Below threshold, the power law re- 
laxation terminates into a pinned configuration, while 
above threshold there is a crossover to linear creep or av- 
erage constant velocity sliding. As it is common for this 
class of systems, the motion of the pileup takes places in 
the form of avalanches whose distribution again can be 
characterized by scaling exponents. 
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II. ELASTICITY 



Two dimensions 



Developing a theory for collective dislocation depin- 
ning requires the basic knowledge of the elastic proper- 
ties of the dislocation assembly in the first place. In this 
section, we determine the elastic response of two par- 
ticular dislocation assemblies: a regularly spaced pileup 
and a low angle grain boundary of edge dislocation lines. 
The two structures are quite similar geometrically; they 
are both one-dimensional arrangements of N dislocation 
lines with the same Burgers vector b and average line 
direction e (for edge dislocations e _L b) along a given di- 
rection of space d, but they differ in the relative orienta- 
tion of the Burgers vector and the arrangement direction 
d. In particular, in a pileup a set of edge dislocations 
lies in the slip plane, defined by the dislocations axis e 
and the Burgers vector, so that d \\ b (see Fig. |21for a 
particular example with e — z and d \\ b \\ y), while 
in the LAGB the edge dislocations lie in the perpendic- 
ular plane such that d J- b (see Fig. O for a particular 
geometry). Neglecting climb, i.e. the motion of a dis- 
location perpendicular to its slip plane, deformations of 
the structure occur solely within the slip plane both for 
the pileup and for the LAGB. In this section we derive 
the shear stress and the elastic energy associated with 
small deformations of these dislocation assemblies. This 
is needed in order to derive the yield stress from statis- 
tical pinning theories. For completeness, we consider the 
problem both in two and in three dimensions. 
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In this particular case, we do not consider elastic de- 
formations along the dislocation line direction. This ap- 
proximation is relevant for thin crystals and amounts to 
treat the dislocations in the structure as rigid lines, ren- 
dering a two-dimensional problem. On the other hand, 
we take into account small variations in the position of 
the dislocations in the one-dimensional structures they 
form. We consider the case of LAGB and then directly 
extend the result to the pileup case. In fact, in linear ap- 
proximation the elastic energy turns out to be the same 
in both cases. 

Here and throughout the paper, we consider an ideal 
LAGB as an infinite set of equally spaced edge disloca- 
tions lying on the yz plane (without loss of generality 
we consider the plane a; = 0) with Burgers vector point- 
ing along the positive x axis h — bx (see Fig. In 
the rigid dislocation approximation, we do not consider 
the deformations along the z axis so that each disloca- 
tion is described by a set of coordinates {xn,yn), where 
yn = nD, D is the grain boundary spacing, and where Xn 
is a small displacement from the x = plane. The shear 
stress field at a givenpoint {x, y) due to a dislocation at 
{xn.yn) is given by jHliil 
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{X - Xn)[ix - Xn)'^ - {V ~ Vnf 



27r(l - V) [{X - Xnf + {y- yn)' 



(1) 

where ^ is the shear modulus and v is the Poisson ratio. 
The glide component of the total force per unit length 
on another dislocation m in the LAGB can be readily 
obtained from the Peach-Koehler expression f = (cr • b) x 
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For small deformations \xm — x„| ^ _D|m — n\ we have 



FIG. 2: A regularly spaced dislocation pileup with Burgers 
vector along the y axis. The ideal configuration is plotted 
with straight dashed lines, whereas the solid lines represent 
their possible glide deformations within the slip plane yz. 
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which can be used to obtain the elastic energy 
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It is instructive to express the elastic energy in Fourier 
space, where one can easily identify the energy cost of the 
different modes. For an infinitely long LAGB TV oo, 
we can write the dislocation displacements as 



BZ 



dk 
2^' 



(5) 



where the integral is restricted to the first Brillouin zone 
(BZ) of the reciprocal space, and using the following re- 
sults 
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we obtain l47 



E ■ 
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(7) 

From this expression, one can see that the elastic in- 
teraction kernel (27r|fc| — Dk^) is not quadratic in the 
wavevector, as it would be the case for a local elastic 
manifold with a constant tension or stiffness, but grows 
roughly as |fc| for long wavelength deformations. This is 
a consequence of long range interactions between dislo- 
cations in the LAGB which render a much stiffer struc- 
ture. In the following sections we will explore the relevant 
consequences of this result in the collective depinning of 
such dislocation structures, something that has been dis- 
regarded in previous dislocation depinning studies. 




FIG. 3: A regularly spaced low angle grain boundary where 
the dislocations Burgers vector is parallel to the x axis. The 
ideal configuration is plotted with straight dashed lines in the 
plane yz, whereas the solid lines represent their possible glide 
deformations within the slip plane xz. 



The elastic energy associated with perturbations of a 
regularly spaced dislocation pileup can be obtained in an 
analogous manner. According to the geometric condi- 
tions assumed here, all Burgers vectors are now oriented 
along the positive y axis, and since the dislocations are 
all in the same slip plane we can now write x = Xn = 0. 
Proceeding as before, the total Peach-Koehler force on 
dislocation m along the new glide direction is given by 



+ 00 



1 



27r(l - u) ^ y„i - yn 



(8) 



Notice that the Peach-Koehler force is now repulsive, 
however the stability of the system is ensured in the case 
of an infinite pileup where the dislocations located at 
the extremes (at ±oo) have fixed positions, or equiva- 
lently, for a finite pileup with periodic boundary con- 
ditions. Thus one can also compute the elastic energy 
cost of small displacements 6ym of the dislocations in the 
pileup with respect to their stable positions. Up to first 
order in 5ym , we obtain a recovery elastic force 
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equivalent to the one obtained for the case of the LAGB. 
The corresponding elastic energy cost 



E 
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(10) 

has the same form with a long-range interaction ker- 
nel. Thus as a partial conclusion of this section, we can 
emphasize that long wavelength distortions of low angle 
grain boundaries and equally spaced pileups of straight 
dislocation lines with translational invariance along the 
dislocation axis have the same nonlocal elastic properties, 
with eigenvalues that grow linearly with the modulus of 
the wavevector considered. 



B. Three dimensions 

In this section, we consider the more general and re- 
alistic case of deformable dislocation lines where, ac- 
cording to the geometry established in the previous sec- 
tion, we loose the translational invariance along the z 
axis. As before, we consider first the case of a LAGB 
with Burgers vectors oriented along the x axis, in which 
each dislocation can be described by a set of coordinates 
(x„(z), y„, z), where again y„ — nD, but where now the 
displacement x„(z) with respect to the a; = plane de- 
pends on the height of the infinitesimal dislocation seg- 
ment considered (see Fig. ^ . The elastic stress field due 
to a general dislocation line or loop can be obtained, for 
instance, by considering the line as being composed of el- 
ementary segments of infinitesimal length jig . Depend- 
ing on the relative orientation of the Burgers vector and 



5 



the local tangent vector f{z), each segment can either 
have edge (f(z) ± b) or screw character (f(z) || b), or 
it can be a combination of both. A first approximation 
of a general dislocation line can be its parametrization 
in terms of a succession of only edge and screw segments 
|48l |. The mathematical form of the elastic stress fields 
generated by these two types of elementary segments is 
simpler and renders amenable the analytic treatment of 
the problem. For instance, the shear stress field due to 
an edge dislocation segment of length Az' and located at 
the point {x',y',z') is given by 46] 



axy{x,y,z) 
where 



fib 



47r(l - v) 
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{y~yr 



Rl 



Az', 
(11) 



Rl = (x-x'f + {v~v'f + {z-z'f. (12) 

On the other hand, the shear stress field due to an screw 
segment of length Ax' is 



, , fib z — z' I 
o-xy{x, y,z) = - — — ^ Aa; . 



(13) 



Equations l(TT|l and ((T^ allow us to calculate the glide 
component of the total Peach-Koehler force f = (a • b) x f 



on an edge or an screw segment. The glide force on an 
edge segment at {xm{z), ym, z) has two components due 
to any other edge f^^ or screw segments f^^ 



fx'^{Xm{z),ym,z) 



fx iXm{z),ym, z) — 
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-^mn {Zj z') 



Az'Az, (14) 



flb^ z — z' dxn{z') 



47r i?^„(z,z') dz' 



Az' Az 



respectively. Note that up to first order in the small 
displacements {xm{z) — Xn{z') — 0), the relative dis- 
tance among segments can be written as R^^{z, z') = 
{Vm — ynY + {z — z'Y . On the other hand, from the gen- 
eral expression for the Peach-Koehler written above, it 
is straightforward to verify that there are no glide forces 
acting upon any screw segment on the dislocation line. 
After summing up all nonvanishing contributions, we can 
obtain the elastic energy as for the two dimensional case 
(see Eq.Q). The elastic energy can be expressed as the 
sum E — E^^ + E^^ of the interaction energies between 
edge-edge and edge-screw segments. These interaction 
energies are given by 
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(15) 
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^mn\Zj Z ) 



(16) 



As we did for the rigid line case, we can also express 
this elastic energy in Fourier space in order to diagonalize 
the interaction matrix and to obtain the wavevector de- 
pendence of the interaction kernel between the different 
deformation modes. The detailed calculation is rather 
lengthy, so we merely indicate the procedure followed and 
the final results obtained. We evaluate separately the en- 
ergy contribution due to the self-interaction between the 
constituent segments of each individual dislocation line, 
i.e. n = m, which we denote by Eq, and the energy con- 
tributions due to the interaction of dislocation segments 



lying on different lines, i.e. n ^ m, which we refer to as 
El . Proceeding this way, we find that the total energy is 
E = E^^ + E^^ + Ef^ + Ef^. We express the disloca- 
tion displacements in terms of their Fourier modes, 

"^"("^ = Lz^l ^e-'^^-e-'?- x(fc, q) (17) 
and evaluate the self-interaction contributions for long 
wavelength deformations qa ^ 1 where a is a short- 
distance cutoff introduced to preclude the interaction of 
a line segment with itself. The result can be written as 
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where 7 is the Euler constant. The usual quadratic 
wavevector dependence of a local interaction kernel is 
modified in this particular case by logarithmic correc- 
tions. This is known to be the result of long range in- 
teractions in the case of a dislocation line, as well as for 



similar singularities such as vortex lines in high temper- 
ature superconductors [ssf . 

Similarly, the energy contributions due to interactions 
between segments of different dislocation lines (n ^ m) 
in the LAGB can be expressed as 
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(21) 



where (^(x) is the Riemann zeta function. Naturally, the 
interaction kernel between the deformation modes for the 
three dimensional grain boundary case depends explicitly 
on both the y and z components of the wavevector in 
an intricate manner. Nevertheless, as in the two dimen- 
sional case, for long wavelength deformations the leading 
term of the interaction kernel is essentially linear in the 
wavevector, which manifests the nonlocality of the inter- 
actions. 

Finally, let us consider the case of a three dimensional 
pileup with the same geometric specifications as adopted 
in the rigid-line case (see Fig. 13). Since the Burgers vec- 
tor is oriented along the y axis, the glide component of 
the Peach-Koehler force on screw segments (now paral- 
lel to the y axis) vanishes. On the other hand, there is 
a non-vanishing glide force on edge segments due to the 
shear stress field generated by any other segments. For 
instance, the shear stress field due to an edge dislocation 
segment of length Az' and located at the point {x' , y' , z') 
in this case is given by ,4^] 



Without loss of generality we consider a pileup ly- 
ing on the X = plane, and we account for small per- 
turbations affecting the ?/„ coordinates of the disloca- 
tion lines in the ideal assembly. Thus, we can replace 
Vm -y-a ^ {Vm + 5ym{z)) " {Vn + hn{z')), whcrc the 
displacements depend on the z coordinate of the infinites- 
imal line segment considered. Expanding up to first order 
in {Sym{z) — 5yn{z')), the resulting Peach-Koehler glide 
forces on an edge segment due to other edge fy^ or screw 
segments fy^ are 



fy^^%y„^ + Sy {z)) = 
lib oy„,[z) - dyn(z ) 



477(1-1/) Rl,n{z,z') 
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(24) 



lib y-y' 



{x~x'f , , 
1-3^^ Az'. 

(22) 

Moreover, the shear stress generated by an screw segment 
of length Ay' is |46| 



where 



Rla,i{z, z') (y™ - ynf + (^ - z'f 
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The corresponding elastic energy is 
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(26) 
(27) 



Proceeding along the same lines indicated in the case of a LAGB, i.e. separating the intra-line interactions 
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from the inter-line ones, it can be expressed in terms of 
Fourier modes as E = E^^ + E^^ + E^^ + Ef^. The 
expressions turn out to be equal to the one obtained for 
the LAGB and can be recovered replacing x{±k, ±q) with 
Sy{±k,±q) in Eqs.UBEIl 

Thus also in this case we find wavevector dependent in- 
teraction kernels whose leading terms (for low wavelength 
deformations) grow either linearly with the wavevector 
in the case of interline interactions, or quadratically with 
logarithmic corrections in the case of intraline interac- 
tions. We thus can conclude that this particular form of 
the elastic kernels is characteristic of the long range inter- 
actions between different dislocations. As we will see in 
the following, the long-range elastic properties preclude 
the analysis of the depinning transition of dislocation as- 
semblies based upon a local theory with a constant effec- 
tive stiffness. 

III. DISORDER: DEPINNING TRANSITION 

Distortions in a LAGB or a pileup arise from interac- 
tions of the dislocations with various kinds of impurities 
such as solute atoms, precipitates or other immobile de- 
fects. The interactions between individual dislocations 
and impurities have been computed and are reported in 
the literature. For the purpose of this article, we will 
consider quenched disorder created by a random distri- 
bution of immobile impurities with concentration c which 
interact with dislocations via a force fp{r) — fogir/^p), 
where /o is the pinning strength, is the interaction 
range and r is the distance between the impurity and the 
dislocation. The detailed shape g{x) of the individual 
pinning force is inessential for most purposes. 

The morphology and dynamics of a pileup or a LAGB 
result from a complicated interplay between elasticity 
and disorder. Pileup and LAGB are examples of the 
general problem of the depinning of elastic manifolds in 
random media, which has been extensively studied in the 
past |14| . In the elastic approximation, the dynamics of 
the pileup or the LAGB follows 

X-Q^ = I d''^x'K{x - x'){u{x') — u{x)) + b(7 + t]{x, u), 

(28) 

where x is a damping constant, a is the applied stress, 
ri{x, u) describes the effect of the pinning centers and the 
elastic interaction kernel K, computed in the previous 
section, scales as in Fourier space. In the following 
we will discuss how the main theoretical approaches to 
the depinning transition can be applied to the problem 
at hand. 



A. Collective pinning theory: weak pinning 

Collective pinning theory describes the behavior of an 
elastic manifold in the limit of weak disorder, when pin- 
ning is due to the fluctuations of the random forces. The 



key concept is the introduction of a characteristic length 
Lc above which pinning becomes effective (or energeti- 
cally advantageous) and consequently the manifold is dis- 
torted. The collective pinning length can be evaluated, 
for instance, by balancing the elastic energy cost and 
the pinning energy gain associated with a small displace- 
ment of a region of linear size L. On scales below L, the 
manifold remains essentially undeformed and, hence, the 
fluctuations in potential energy follow Poissonian statis- 
tics. The effective concentration of the pinning defects 
along the LAGB is given by 

The first expression refers to pinning by columnar defects 
of areal concentration c in d = 2, and the second to 
pinning by localized defects of volume concentration c in 
d — i. In d=2, the characteristic energy of a section of a 
LAGB of size L displaced by an amount of the order of 
u can be written as 

J 2 2 

^ ^ - Mp^c^sLu. (30) 

Here both E and /o are defined as quantities per unit 
length. In the case of a thin film of thickness /i, one 
can obtain their three-dimensional counterparts just as 
E — hE and /o — hfo. Note the scale-independence of 
the nonlocal expression of the elastic energy fj,b^u^/D^ 
in contrast to what would be this energy in the local ap- 
proximation fib^u^ / DL. Essentially the same expression 
holds for the pileup. Balancing elastic and pinning con- 
tributions and imposing that the displacement is of the 
order of the pinning range w ~ ^p, one readily obtains 
= {fi^b^^p)/{D^f§Cee). The LAGB is depinned when 
the work done by the external stress in moving a segment 
of length Lc over the distance exceeds the character- 
istic pinning energy E{Lc) of this segment. Equating 
E{Lc) = acbLc£,p/ D, for the case above the result is given 
by aab = {c,sf^D^)/{fib^). 

A similar calculation in d = 3 is more subtle, since 
the elastic and the pinning energies scale with the same 
power of L and thus cancel in the simple dimensional 
approach discussed above. As we will discuss in the next 
section, this reflects the fact that d = 3 is the upper 
critical dimension for the transition. To obtain Lc in this 
case, one should perform a perturbation expansion in the 
disorder, as discussed in Ref. |0| in the context of the 
flux line lattice. One essentially computes the typical 
displacement u for a system of size |r| = L, which for a 
LAGB is given by 



(|.M-..(0)R^/|^/|i;(i-».k,) 

G(k)G(k')F(k)i^(k') (31) 
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where G(k) is the Green function associated to the elas- 
tic kernel determined in the previous section, and i^(k) 
is the pinning force density. In the spirit of collec- 
tive pinni ng the ory (F(k)i^(k')) = WS'-^'^ {k + k') with 
W = (/o\/coffCp)^- The explicit calculation leads to the 
characteristic displacement 



(32) 



This expression can then be inverted, imposing u ~ ^p, 
to obtain 



Lc~D exp 



(33) 



The depinning stress can then be obtained as in d = 2 
and is given by acb = {fj,b^S,p)/{DLc)- Again these re- 
sults generalize directly to the pileup case. It is however 
important to note that they refer to the continuum limit, 
when one can neglect the discrete nature of the disloca- 
tion system. To be consistent with this assumption one 
should have D. 



B. Strong pinning: Friedel statistics 

Collective pinning is due to a statistical superposition 
of the forces created by many obstacles. In the limit 
of strong and/or diluted pinning centers, however, the 
characteristic bulge of width and extension Lc as en- 
visaged in the previous section may not interact with 
enough pinning centres for this viewpoint to be valid. 
Simple estimates for the boundaries of the collective pin- 
ning regime are given by the inequalities LcCp ^ l/^off 
and L^^p > l/coS for the d=2 and d=3 cases discussed 
above, respectively. 

In the regime of strong pinning, dislocations are pinned 
by individual obstacles. The spacing of obstacles along 
the dislocation and the depinning stress can be obtained 
by an argument which was, in the context of single dis- 
locations, developed by Friedel. The basic idea is to con- 
sider the behavior of a dislocation segment as it depins 
from a pair of strong obstacles. The length of the seg- 
ment is L, and it forms a bulge of width u. If the dis- 
location segment overcomes one of the pins it will travel 
by an amount which is, again, of the order of u and, 
hence, sweep an area of the order of Lu. Now we can 
estimate the depinning threshold by requiring that dur- 
ing this process the freed dislocation segment encoun- 
ters, on average, precisely one new obstacle. In other 
words, precisely at the point of depinning the disloca- 
tion starts to move through a sequence of statistically 
equivalent configurations. For a dislocation this leads 
to the condition Lu ~ l/(c^p). L and u can be re- 
lated by equating the work done by the external stress 
a in bulging out the dislocation to the concomitant elas- 
tic energy increase, Vv? / L = abuL, where F is a con- 
stant line tension. Finally, the depinning force can be 



obtained by comparing the external force haL with the 
pinning force /q. Solving these three equations, one ob- 
tains the Friedel length -Z^/ — (r/c^p/o)"^^^ and the de- 
pinning stress Gch ~ (c^p/p/F)^/^. 

This argument can be generalized in a straightforward 
manner to the case of grain boundaries. Let us first con- 
sider the depinning of a two dimensional LAGB as dis- 
cussed above in the weak pinning limit: In this case, the 
Friedel condition reads Lu ~ 1/coff, the elastic energy 
per unit length of a bulge of width u and extension L 
is ^b^u^ / which must equal the work per unit length 
abLu/D, and the force balance (again per unit length) 
is acbL/D = /q. Combining these relations we find that 
the Friedel length and the depinning stress are 



L 



f 



^ib^ 



nb^ 



(2d). 



(34) 



In 3d the Friedel condition is L^u ~ l/ccs, the energy 
balance reads nb^u^L/D^ = abL^u/D, and the force bal- 
ance is acbL'^/D = /q. This yields 



L 



f 



fib^ 



D'^Ccsfa 



„2 f3r)5 



(3d). 



(35) 



Table 1 presents a compilation of results for the weak 
and strong pinning cases in two and three dimensions. 
For comparison we have also included results obtained 
under the assumption that the elastic behavior of the 
grain boundary can, in local elasticity approximation, be 
described by a scale-independent surface energy Fq ~ 
^ib^/D. 



C. Zener pinning in grain growth 

^From an experimental viewpoint, grain boundary pin- 
ning is important since the mobility of grain boundaries 
may be a limiting factor in grain growth [40l l4ll | . The 
dependence of the average grain size R on the impurity 
concentration c goes under the name of Zener formula 
and was originally proposed in Ref. |4Clj | and generalized 
afterwards j41| 



R cc c 



(36) 



where m is the Zener exponent. Several theoretical and 
numerical derivations of this formula have been discussed 
in the literature, often based on the local elasticity ap- 
proximation '41', '42^ . 

Grain growth is driven by a reduction in energy: For an 
average grain size R and straight grain boundaries, the 
characteristic energy stored per unit volume in the form 
of GB dislocations is of the order of Fo/i? and, hence, the 
energy gain achieved by increasing the grain size by di? is 
(Fo/i?^) di?. Physically, the removal of GB dislocations 
occurs through the motion of junction points in the GB 
network. As junction points must drag the connecting 
GB with them, which may be pinned by disorder, motion 
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model dimension 


type of elasticity 


type of pinning 


pinnin 


g length 


critical stress 


Zener exponent 


2d 


local 


weak 






1/3 






m = 2/3 






2d 


local 


strong 


Lf = 


J- 


1/2 


o-cb = 




j 1/2 


m = 1/2 


V /O^cff j 


I To 


2d 


non-local 


weak 


= - 






(7c6 = 






m — 1 




To 


2d 


non-local 


strong 


Lf = 


To 




(7c6 = 






m — 1 




To 


3d 


local 


weak 




( rj^Sp ^ 


1/i 




m — 1 


3d 


local 


strong 


Lf = 


( ro > 

\ /oCcff J 


1/2 




m — 1 


3d 


non-local 


weak 


Lc — D exp 




(JcO = 


^exp 




exponential 


3d 


non-local 


strong 


Lf = 


To 




n 


m = 2 



TABLE L Overview of pinning stresses and pinning lengths obtained from different models 



can only occur if the energy gain at least matches the 
dissipative work which has to be done against the pinning 
forces. The dissipative work per unit volume expended 
in moving all GB by dR is aj)/ {DR)AR^ and balancing 
against the energy gain yields the limit grain size 

(37) 

According to this relation, the grain size is inversely pro- 
portional to the pinning stress. This gives a possibility to 
experimentally assess the nature of the pinning by 'tun- 
ing' the pinning stress and measuring the grain size as a 
function of the tuning parameters. 

An obvious method to tune the pinning stress is to 
modify the concentration of the pinning centers and to 
measure the impact this has on grain size 40]. However, 
this requires the comparison of results from different sam- 
ples and experimental results reported in the literature 
are inconclusive 49, 50] . We therefore refer to a quite 
different and unusual type of grain growth experiment 
where the configuration of the pinning obstacles is kept 
constant but the properties of the lattice are changed. 

This type of grain growth experiment may be carried 
out on vortex lattices of type-II superconductors in which 
quasi two-dimensional grain structures are observed in 
(vortex polycrystals). In such a system an external mag- 
netic field penetrating the sample forms flux lines that 
are disposed in a triangular lattice, whose elastic proper- 
ties (namely the lattice spacing and the elastic moduli) 
depend on the magnetic field itself. A detailed theory 
of grain growth in vortex polycrystals [s^ ] can be devel- 
oped along the same lines followed here. In this case, 
it is worth pointing out that the agreement of the the- 
ory with magnetic decoration data for the average grain 
size is quite satisfactory, especially if compared to the 
estimate based on local elasticity assumptions. 



IV. DYNAMICS: CRITICAL SCALING NEAR 
THE DEPINNING THRESHOLD 

Second-order phase transitions can be described by 
scaling laws and critical exponents and the depinning 
transition is no exception. In the system discussed here, 
the control parameter is the applied stress, so that scal- 
ing laws depend on the distance from the critical point 
a — Uc- In particular, as the system approaches the tran- 
sition the correlation length diverges as ^ ~ (ct — ac)~'^ ■ 
Similarly, one can define a characteristic correlation time 
t* , related to the correlation length as t* ~ The 
average dislocation velocity reaches a steady value, scal- 
ing as V {a — (Jc)^ , above the transition, and vanishes 
below. Before the steady-state the average velocity de- 
cays as a power law t^^ , for times t < t* . Furthermore, 
the Orowan relation, a phenomenological law which re- 
lates the rate of plastic deformation 7 to average values 
of moving dislocations density and velocity in a given 
crystal, implies that similar scaling laws should hold for 
the strain rate 7 = bpv, where p is the density of mov- 
ing dislocations, and v their average velocity. In this 
respect, it is thus tempting to establish a relationship be- 
tween dynamical behavior and the creep laws observed in 
plastically deforming crystals, i.e. the crossover between 
primary (power law) to secondary (linear) creep, due to 
their resemblance. 

Scaling exponents also characterize the morphology of 
the dislocation arrangement, which exhibits roughening 
close to the depinning transition. The roughness can 
be quantified measuring the average displacement cor- 
relations C{x — x') = {{u{x) — u{x')Y). At the transi- 
tion in the steady-state, we expect a self-affine scaling 
C{x) ~ x^^, where C, is the roughness exponent while 
the transient behavior is described by a scaling form of 
the type C{x,t) = t^* f{x/t^^^). As in ordinary criti- 
cal phenomena, only a fraction of the scaling exponents 
are independent. For instance, one can easily derive the 
relations Pt — z/C, and 9 — (3/{vz). 

In general, it has been shown that in the depinning 
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problem there are only two independent exponents that 
have been computed using the renormalization group. To 
connect our problem to previously obtained results, we 
notice that the effective elastic energy of the pileup and 
LAGB scales as \q\ in Fourier space, as in the problems 
of contact line 51] and planar crack depinning |^. We 
can thus directly apply to our case the results obtained 
for a manifold with long-range elastic energy [ij, l2C|. 
The renormalization group analysis predicts that dc = 5 
is the upper critical dimension, above which fluctuations 
are suppressed. Thus for d > there is no roughening 
(i.e. C = 0) and the other exponents can be computed in 
the mean- field approximation, yielding (3 = z = v = 1. 
These results are valid in the physically interesting di- 
mension d — 3 apart from additional logarithmic correc- 
tions. For d < 3, a renormalization group expansion in 
e — 3 ~ d has been performed to compute the exponents 
which at first order in e are given by f3 = 7/9, v = 3/2, 
( — 1/3 and z = 7/9 Using the scaHng relation 

9 — j3/{vz) one obtains = 2/3, which coincides with the 
exponent of the so-called Andrade creep law, observed in 
the creep deformation of several materials [1^ |43| ■ 



A. Two-dimensional pileup: Numerical simulations 

In order to test the validity of the theoretical results, 
we have performed a series of numerical simulations for 
the dynamics of a two-dimensional pileup. This corre- 
sponds to an effective one dimensional model in which N 
interacting point dislocations move along a line in pres- 
ence of quenched disorder. For simplicity we consider 
periodic boundary conditions, so that in absence of dis- 
order the equilibrium configuration is an equally spaced 
pileup. To test the dependence on the system size, we 
change the dislocation number TV and the system size L 
keeping the dislocation spacing D = L/N constant. 
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FIG. 4: The decay of the average pileup velocities as a func- 
tion of the applied stress a. For a > — 0.675 the velocity 
reaches a steady value and decays to zero otherwise. 



The equation of motion for the dislocation i in the 
pileup is given by 



dxi 
dt 



Aifo^^^^— +6<7 + ^/p(a;,-Xp), (38) 



where x is an effective viscosity and a is the applied 
stress. The pinning centers are placed at randomly cho- 
sen positions Xp (with P — l,...Np) and exert an at- 
tractive force on the dislocations 



-Jo^e 



(39) 



In order to correctly take into account the effect of peri- 
odic boundary conditions the interactions between dislo- 
cations are summed over the images. In one dimension 
the sum can be performed exactly and l/|x| in Eg. 1381 is 
replaced by 



E 

k—~oo 



1 



X + kL Lt'An{Tix/L) 



(40) 



The equation of motion (Eq. I38|l is integrated numeri- 
cally using a Runge-Kutta algorithm for different values 
of the applied stress. We take as initial condition a per- 
fectly ordered pileup, with equally spaced dislocations. 
For the simulations reported here, we first considered 
N = 64, 128, 256, 512 dislocations with a spacing D — IQ 
and average pinning center spacing dp = L/Np = 2. 
The units of time, space, and forces are chosen so that 
fib^ = 1, X = 1, and 6=1, and we set /o = 1 and = 1. 
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FIG. 5: The decay of the velocity aA^ a > Oc — 0.675 for 
different values of A*'. As N increases the power law scaling 
region extends. The line has a slope of 6^ = 0.65. 

In Fig. ^ we report the time decay of the average pileup 
velocity for different values of the applied stress. For 
large stress values, a > Uc — 0.675, the initial power law 
decay is followed by a plateau, while the velocity decays 
to zero otherwise. This allows to identify the depinning 
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point as <Jc — 0.675. This is confirmed by the finite size 
analysis shown in Fig. |S1 indicating that for Uc = 0.675 
the power law extends further as the system size is in- 
creased. The exponent of the power law scaling 6 ~ 0.65 
is in good agreement with the theoretical expectations. 

Moreover, in order to characterize the growth of corre- 
lations at the critical point, we compute the displacement 
correlation function C{i — j,t) = {{{ui{t) — Uj{t))'^)y/^ 
at different times i for cr = CTc (see Fig. The 
curves can be collapsed using the scaling form C{x,t) = 
t^/^/(a:/ti/^) with C = 0.35 and z = 0.9 (see the inset 
of Fig. IHJ). To confirm this result we have also com- 
puted the evolution of the power spectrum P{k,t) = 
J dxC{x)ex.p{ikx). These curves can also be coUapsed 
as P{k,t) = ^'^^^"'"^'/^(/(fci^/^) with the same exponent 
values as the correlation function. 




1 10 100 1000 

D(i-j) 

FIG. 6: The growth of the correlation function at the depin- 
ning transition at different times. The data collapse in the 
inset allows to estimate the roughness exponent — 0.35 and 
the dynamic exponent z = 0.9. 

In summary, all the exponents determined from the 
simulations are in good agreement with the renormal- 
ization group predictions and with previous simulations 
based directly on the elastic approximation, confirming 
the validity of the elastic theory for the pileup. 

B. Three-dimensional pileup: relaxation of 
slip-band growth rates 

Direct experimental observation of the velocity relax- 
ation behavior of planar dislocation arrangements is pos- 
sible in certain alloys exhibiting so-called planar slip 
where dislocations may form huge pile-ups (see Figure^). 
The motion of these planar dislocation groups goes along 
with the formation of large slip steps along the traces 
where the slip plane of the pileup intersects the surface 
of the metal specimen. For a moving pile- up consisting of 
roughly equally spaced dislocations, the slip step growth 




kt 



FIG. 7: The power spectrum of the pileup at the depinning 
transition. The data collapse is consistent with the scaling of 
the correlation function. 



rate is proportional to the dislocation velocity. Since the 
slip step growth rate can be directly observed, this gives 
a possibility for a direct experimental check of theoretical 
predictions. 




time [s] 

FIG. 8: Growth rate of slip steps on the surface of Cu-30at% 
Zn deformed at room temperature as a function of the time 
passed after growth has started; after Ref. 13- The line is a 
power law with exponent 9 = 1. 

Figure |H1 shows observed rates of slip step growth 
as a function of the time after growth has started. 
The double-logarithmic plot indicates relaxation of the 
growth rate (the dislocation velocity) according to w oc 
with a characteristic exponent 6 = 1 ± 0.1 over six 
decades. This indicates that the power-law is related to 
a depinning transition of planar dislocation arrays in 3D 
for which we expect mean-field exponents (see above) 
and according to the scaling relation 9 = (3/{vz) a value 
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— 1. The apparent length of the scaUng regime in- 
dicates that driving of the dislocation arrays occurs at 
stresses very close to the critical one. This is in line with 
the general observation that dislocation arrangements in 
slowly deforming crystals (where 'slow' covers the entire 
range of strain rates used in typical experiments, 01) 
are in a close-to-critical state ISJ, ISal . 



V. CONCLUSIONS 

We have investigated the depinning transition of pla- 
nar dislocation arrays such as small-angle grain bound- 
aries or dislocation pile-ups. Contrary to the case of iso- 
lated dislocations, the elastic interactions between dis- 
location line segments in such arrays are of long-range 
nature and, hence, cannot be described within a line- 
or surface-tension approximation. The pinning of planar 
dislocation arrays has been investigated both in the weak 
and strong pinning limits using collective pinning theory 
and Friedel statistics, respectively. The results have been 
applied to the problem of grain growth limited by grain 
boundary pinning referred to as Zener pinning. 

Long-range elastic interactions also govern the dynam- 
ics of planar dislocation arrays at the depinning thresh- 
old. In two dimensions, computer simulations and the- 
oretical arguments suggest that the dynamics falls into 
the same class as contact-line depinning, while in three 
dimensions the dynamical behavior can be described by 
mean-field exponents. In particular, we have demon- 



strated that the mean-field prediction for the veloc- 
ity relaxation of a planar dislocation array (pile-up) is 
consistent with experimental observations of the time- 
dependent growth of slip bands in alloys exhibiting pla- 
nar slip. 

The dislocation arrangements discussed in the present 
study have a simple, quasi-planar geometry in which only 
dislocations of one sign are present and only small pertur- 
bations of the planar arrangement of the dislocations are 
permitted. Because of this particular geometry, the dis- 
location assemblies behave like two or three- dimensional 
long-range elastic manifolds. The situation is much more 
complicated when dislocations of different types and di- 
rections of motion have to be considered. In such sit- 
uations, there is still a transition between a stationary 
and a moving state of the dislocation assembly ('yielding 
transition'). However, in general dislocation assemblies 
the existence of metastable stationary states does not 
depend on the presence of quenched disorder as in the 
present study. Rather, the interactions between dislo- 
cation lines of different type together with the dynam- 
ics constraints which tie the motion of the dislocation 
lines to their respective slip planes lead to the possibil- 
ity of forming metastable jammed configurations even in 
the absence of any disorder. While the general scenario 
of dynamic non-equilibrium phase transitions applies to 
such systems, no ready-made theoretical framework is 
available and, hence, a theory of the yielding and dy- 
namic behavior of general dislocation systems remains a 
formidable task for future investigations. 
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